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Because the physical interpretation of the spin-polarization of a ferromagnet determined by point- 
contact Andreev reflection (PCAR) is non-trivial, we have carried out parameter-free calculations 
Q ' of PCAR spectra based upon a scattering-theory formulation of Andreev reflection generalized to 

CZ3 . spin-polarized systems and a tight-binding linear muffin tin orbital method for calculating the corre- 

sponding scattering matrices. PCAR is found to measure the spin-dependent interface transparency 
rather than the bulk polarization of the ferromagnet which is strongly overestimated by free electron 
model fitting. 
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The spin-polarization of a ferromagnet at the Fermi energy (not to be confused with the total magnetic mo- 
O , ment) is a parameter of crucial importance entering all models of magneto-electronic phenomena such as tunneling 
magnetoresistance (TMR). Although Moodera et aL's TMR experimentsa could be interpreted in terms of bulk spin- 
y—{ \ polarizations using a very simple model due tOpJullierea, subsequent experiments on a wider variety of materials 
■ cast doubt on the general validity of the modelfTTl They also led to the need to find simpler ways of measuring 
" spin-polarizations than the benchmark determinations by Tedrow and MerseveyQ using tunneling spectroscopy at 
ferromagnet | insulator | superconductor (S|I|F) junctions. These are difficult to make and characterize, especially for 
complex materials. Andreev reflection spectroscopy at ferromagnet-superconductor (F|S) interfaces in ballistic point 
contacts (PCAR) has been proposed as a simpler and complementary method to measure spin polarizationsQcl. When 
^vq ■ the diameter of the point contact is smaller than the mean free path but larger than the Fermi wavelength, transport 
is governed by the Sharvin resistanceEJ, appropriately modified to take into account the interface "transparency" . 
In the Blonder-Tinkham-Klapwijk (BTK) theoryM, this is done by representing the interface as a planar <5-function 
barrier between free-electron materials. The Bogoliubov-de Gennes equationEj for a normal-superconducting (N|S) 
interface is readily solved for this model and the essential physics expressed in terms of a dimensionless barrier strength 
Z and the superconducting energy gap A. At the interface to a ferromagnetic material, the Andreev reflection is 
suppressed by the difference between majority and minority electronic band structurest 3 ! Using the parabolic-band 
Stoner-Wohlfarth model for ferromagnetic materials, BTK theory can easily be generalized to spin-polarized systems. 
\ With only a single additional parameter P, the spin polarization of the ferromagnet, a good fit to experimental data 
r* ■ is apparently achievedoEI. 
, | In spite of this, doubts have been raised about the definition and-physical significance of the parameter P derived 
from different experiments for all but 100% polarized ferromagnetalj. Apart from concerns about the relevant Fermi 
surface averaged quantity to be used to describe different transport regimes, transport through an interface should not 
depend solely on parameters characteristic of bulk materials. The representation of an interface between materials with 
quite different electronic structures as a (5-function barrier can only be justified in terms of analytical simplicity. The 
free-electron Stoner model does not describe itinerant ferromagnetism in transition metals in an internally consistent 
manner and the probability of Andreev reflection is overestimated by such models; complex transition metal Fermi 
surfaces lead to a reduction of the phase space available to an electron to pair with a . sp in-flipped retro-reflected hole. 

In this letter we present a generalized scattering formalism of Andreev reflectiont3ll3 to spin-polarized systems and 
calculate the PCAR spectra with and without interface disorder. By basing our study on ab-initio calculations we 
can treat the full complexity of the transition metal band structure without introducing any arbitrary parameters 
and make contact with experiments for a number of systems of current experimental interest. We will be forced to 
conclude that the generalized BTK model omits an important aspect of the problem, basically because the bimodal 
distribution of the transmission probabilities-.of real interfaces are inaccurately modelled using a (5-function potential. 
We start with the Bogoliubov-de GennesEj equation for the electron and hole wave functions of a superconductor 
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where a = ±1. Hq^ is the single-electron Hamiltonian matrix for majority (<r = 1) and minority (er = —1) spins. 
Cecr(h-o-) are the coefhcient-vectors of the wave functions of electrons (holes) in some convenient basis. The excitation 
energy e is measured relative to the Fermi energy. Following BeenakkeitJ, wave-function matching at the F|S interface 
is achieved by inserting between F and S a fictitious region in which S assumes its normal state (disregarding the 
proximity effect). At the F|N interface, the scattering of the Bloch states for electrons and holes can be written 



/ C~(F) \ 

ctAN) 

C h-a( F ) 



S a (e) 
S- a *{-e) 



( C+(F) \ 



(2) 



where +(— ) denotes right (left) going waves and 



rfi(e) *? 2 (e) 
^i(e) »fa(e) 



(3) 



is the normal-state scattering matrix. The subscript 1 refers to F, 2 to N. At the (fictitious) N|S interface there is 
only Andreev scattering: 

C- a {N) = aC^_ <J {N)e i <t>, (4) 

C+_ a (N) = a*C+(N)e-^, (5) 

where for |e| < Ao, a = exp [— i arccos(e/Ao)] describes the phase shift due to the penetration of the wave function 
into the superconductor and for \e\ > Aq, a — (e — sgn (e) y/e 2 — Aq) /Aq, \a\ < 1. The incoming and reflected waves 
are related on the F side by 
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are the reflection coefficients for the total system. In terms of the spectral conductance 
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and the Fermi-Dirac distribution function /, the current-voltage relation is eI(V) = J de[f(e) — f(e + eV)]G fs{ £ )- 
The charge transport in the presence of Andreev reflection is thus expressed in terms of th£ scattering matrix for the 
normal state. This expression is the spin polarized generalization of Beenakker's formulalij and reduces to the BTK 
equivalent if we use the scattering matrix of a <5-function potential and the parabolic band model. Note that the 
expressions also hold for energies above the superconducting gap. 

The parametep-free calculation of the reflection matriceal3 and the conductance (|9)-js based on the surface Green's 
function methodic implemented with a tight-binding linear muffin tin orbital basistH Because a minimal basis set 
is used, we are able to carry out calculations for large lateral supercells and model disorder very flexibly within 
such supercells without using any adjustable parameters. The electronic structure is determined self-consistently 
within the loc&Uspin density approximation. For disordered layers the potentials are calculated using the layer CPA 
approximatior£3. n 

We focus our attention here on the Cu/Pb, Co/Pb and Ni/Pb systems investigated experimentally inQ which were 
grown with an /cc(lll) orientation in point contacts. We encounter the practical difficulty that the lattice constants 
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of Pb are quite different from those of Cu, Co and Ni. Lattice matching is modelled using lateral supercells containing 
4x4 Cu (Co or Ni) atoms on one side of the interface and 3x3 Pb atoms/monolayer on the other side or 7 x 7 
Cu, Co or Ni and 5x5 Pb atoms/monolayer which gives better matching but is computationally more expensive. 
Disorder is taken into account by introducing a two monolayer-thick 50% interface alloycJ. The results are not very 
sensitive to moderate variations in the alloy concentration in this range. Experimental! have been carried out at 
slightly elevated temperatures to suppress a possible proximity effect. Our computations are therefore carried out 
using the corresponding values of kT in the Fermi function: 4.2 K for Cu/Pb, 1.4 K for Co/Pb and 2.5 K for Ni/Pb. 
The only parameter in the calculations is the superconducting energy gap Apb which is taken from experiment and 
is very close to the bulk value. 

For an ideal N|S interface Andreev reflection leads to a doubling of the conductance. In Table | the zero bias 
values Gfs(0)/Gfn{0), calculated for clean and alloy interfaces, are compared with the experimental values of 
Rcf.El. Suppression of the Andreev reflection is apparent in all three cases but is smallest for non-ferromagnetic 
Cu. Even there, however, bandstructure mismatch between Cu and Pb gives rise to a significant reduction of the 
interface transparency and suppression of AR. The agreement with experiment is remarkably good and only in the 
case of Co/Pb is there any indication of a discrepancy. This picture changes when we compare the calculated and 
experimental spectra in Fig. |l| shown in terms of g(V) = [Gfs{V) — Gfn{0)] /Gfn{0) where Gfn(0) is the normal 
state conductance and Gfs(V) is the differential conductance for an F|S interface at finite temperature and finite 
bias. 

Let us first focus on the paramagnetic Cu/Pb point contact, Fig. |l|(a). Whereas results for a specular interface 
deviate significantly from the measured data, the decreased conductance of a rough interface leads to essentially 
perfect agreement with experiment. Aspects of the problem not taken into account, like spin-orbit scattering and 
strong-coupling superconductivity in Pb are apparently not very important. 

The agreement with the PCAR spectra for the ferromagnetic systems (Figs. |l](b) and (c)) is not satisfactory close 
to the band edge of quasi-particle excitations, where the theoretical results appear to strongly exaggerate the Andreev 
reflectionEl Since we have no parameters with which we could improve the agreement the deviations should be sought 
in deficiencies of the model. For example, the closa.proximity of the ferromagnet might smear out the peaked density 
of states of the superconductor at the band edgescil. Whereas this point deserves attention, we focus in the following 
on the zero-bias conductance. 

The spin-polarization of the interface conductance, P = (G rna j — G m i n )/ {G ma j + G m i n ), will in general differ from 
that of the bulk ferromagnetic material. For example, for Ni/Pb, P = 0.5%,-for a specular interface and —1.4% for 
a rough interface. This should be compared to the BTK fit obtained in RefE of Pbtk — 32%. Our prediction of a 
small spin-dependence of the interface conductance could be tested directly by measuring the magnetoresistance of 
a Ni/Pb/Ni trilayer (or multilayer) in the current-perpendicular-to-the-plane configuration with the layer thickness 
chosen sufficiently thin that the total resistance is dominated by the interface contribution. Using the expression given 
by Schep et a/£3 for diffuse systems, we predict interface resistances of R ma j — 0.97 fQm 2 and R m in — 1-27 fflm 2 
for a rough Ni/Pb (111) interface. 

We emphasize the difference between the single-channel BTK model and the multichannel scattering theory gener- 
alization by focusing on the zero-voltage, zero-temperature conductance (^) for a non-magnetic system: 
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is expressed in terms of the eigenvalues of the transmission matrix, T„, and the sum is over the propagating channelsEi. 
For the Cu/Pb interface, the distribution function of transmission matrix eigenvalues, p(T), is shown in Fig. 0(d). It 
turns out to be a bi-modal function with a clear dependence on interface roughness and on the constituent materials. 
Obviously, p(T) strongly affects the Andreev reflection probability with the main contribution coming from the highly 
transmitting part, whereas the average conductance is much more affected by the low transmission channels. Using 
the Andreev conductance to derive an average transmission is therefore misleading. For a ^-function potential and free 

electrons the transmission probability as a function of the in-plane wave vector has the form = 1/ ^1 + i-(^/k F ) 2 ) ■ 

The distribution function for Z = 0.3, the value which must be taken to reproduce the average conductance, is plotted 
as a dashed curve in Fig. 0(d) and is seen to be qualitatively incorrect compared to the realistic distributions. 

For the F|S interface the situation is more complicated. Since Cooper pairs are phase sensitive, the conductance is 
governed by the full scattering matrix, i.e. by the scattering amplitudes and not just the transmission probabilities. 
We therefore cannot use a single distribution function to characterize the interface because the Andreev reflection 
probability contains terms like 1 — Q^r^i^T^* (~ e ) which cannot be expressed by Tj and Tj. Still, it should be clear 
that the deficiencies of the 6 function barrier potential will not vanish. 

In summary we calculate the PCAR spectra of Co/Pb, Ni/Pb, and Cu/Pb point contacts and the spin-polarization 
of Co/Pb, Ni/Pb interfaces from first principles which allows us to remove what is a dubious approximation made 
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when applying the BTK model to transition metals. We find good agreement with experimental results by the Cornell 
group for the full spectrum of Cu/Pb and for the zero bias conductance in the ferromagnetic systems. Qualitative 
differences are revealed by comparing our results with that of a simple BTK model. The physical significance of 
parameters obtained by fitting experiments with this simple model is doubtful. 
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Clean Interface 


1.54 


1.29 


1.08 


Alloy Interface 


1.36 


1.16 


1.00 


Experiment 


1.38 


1.18 


1.13 



TABLE I. Zero bias suppression of Andreev reflection compared with the experimental values obtained from Ref. |^. 
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FIG. 1. (a)- (c) Calculated g(V) curves for disordered (solid line) and clean interfaces (dashed line) compared with the Cornell 
group's experimental results (o). (a) Cu-Pb, (b) Co-Pb, (c) Ni-Pb. (d) Distribution function of the transmission matrix 
eigenvalues of a clean (x) and disordered (o) Cu/Pb interface. The distribution function of the BTK model with Z = 0.3 is 
plotted as a dashed line for comparision. 
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